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Searches for dark matter production at particle colliders are complementary to direct-detection 
and indirect-detection experiments, and especially powerful for small masses, < 100 GeV. An 
important collider dark matter signature is due to the production of a pair of these invisible par- 
ticles with the initial-state radiation of a standard model particle. Currently, collider searches use 
individual and nearly orthogonal final states to search for initial-state jets, photons or massive gauge 
bosons. We combine these results across final states and across experiments to give the strongest 
current collider- based limits in the context of effective field theories, and map these to limits on 
dark matter interactions with nuclei and to dark matter self-annhiliation. 



PACS numbers: 

Though the presence of dark matter in the universe 22 
has been well-established, little is known of its particle 23 
nature or its non-gravitational interactions. A vibrant 24 
experimental program is searching for a weakly interact- 25 
ing massive particle (WIMP), denoted as Xi and inter- 26 
actions with standard model particles via some as-yet- 27 
unknown mediator. If the mediator is too heavy to be 28 
resolved, the interaction can be modeled as an effective 29 
field theory with a four-point interaction. 30 

One critical component of this program is the search 
for pair-production of WIMPs at particle colliders, specif- 
ically pp — > XX the LHC via some unknown interme- ^3 
diate state. As the final state WIMPs are invisible to 
the detectors, the events can only be seen if there is as- 
sociated initial-state radiation of a standard model par- 
tide flMI , see Fig [I] recoiling against the dark matter 37 
pair. 38 

The LHC collaborations have reported limits on the 
cross section of pp — >■ XX + where X is a gluon or 
quark [H [5] , photon [HJ [7] , and other searches have been " 
repurposed to study the cases where X is a [8] or Z 

43 
44 
45 

g, 7, Z, or W 




53 

FIG. 1; Pair production of WIMPs (xx) in proton-proton 54 
collisions at the LHC via an unknown intermediate state, with 
initial-state radiation of a standard model particle. 



boson [9l|T0]. In each case, limits are reported in terms of 
the mass scale of the unknown interaction expressed 
in an effective field theory [THSl [TIHT7] . These various 
initial-state tags probe the same effective theory, but are 
largely statistically independent due to their nearly or- 
thogonal event selection requirements. As the relative 
rates of radiation of gluons (quarks), photons, W oi Z 
bosons from the incoming quark (gluon) legs are deter- 
mined by the standard model, the various probes may be 
combined to give the strongest limits without any loss of 
generality or additional theoretical assumptions. 

Recently, an analysis of multi-jet final states was shown 
to add some sensitivity to the mono-jet analyses [TF|; that 
sample is not statistically independent from the mono-jet 
results used here, and is not included. An earlier global 
analysis of indirect and direct constraints with Tevatron 
data and mono-jet data from ATLAS provided an initial 
set of combined constraints [19, using the approximations 
of a x^ technique. 

In this paper, we perform a full statistical combina- 
tion of the limits from all available channels (mono-jet, 
mono-photon, mono-Z^) from both ATLAS and CMS at 
^/s = 7 TeV, accounting for the dominant correlations 
and providing the most powerful current collider con- 
straints. While the limits reported by the experimental 
collaborations are typically given for a few select effec- 
tive operators, we calculate the efficiencies of their selec- 
tions and reinterpret their searches for the complete set 
of operators relevant for Dirac fermion or complex scalar 
WIMPs. 

Models 

The effective theories of dark matter considered here 
consider the possibility that the final-state WIMPs are a 
Dirac fermion (operators D1-D14 in Ref [T3]) or a com- 



Final states with a heavy boson have little power relative to 
mono-photon or mono-jet; we include mono-Z as a demonstra- 
tion, and do not include mono-VF, although see [8]. 
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FIG. 2: Cross sections for pp — > XX + ^ production where X 
is initial-state radiation of a jet, photon or Z boson. Jet and 
photon final states include a pT > 80 GeV cut at parton-level. ' 
Each pane shows the cross section for a different effective 
operator: top is D5, center is D8, bottom is D9. See Ref. [3] c 
for operator definitions. , 



plex scalar (operators C1-C6 in Ref [Ml). These four-™ 
point effective operators assume that the unknown in-^"^ 
termediate particles have a heavy mass scale; we use a 
suppression scale, M^,. Cross sections at leading order 
for production in pp collisions at ^/s = 7 TeV are shown 
in Fig [2] for select operators with — 1 TeV for il-„4 
lustration. Recently, next-to-leading-order calculationSios 
have been performed for mono-jet and mono-photon pro- 
cesses (201 showing ratios of crNLo/fLO ~ 1-2 — 1.5; ourme 
mono-jet results partially include this effect by generat-w? 
ing and matching multiple-parton emission. 

108 

For some operators, cross sections of dark matter pro- 
duction at the LHC can be transformed into cross sec-ios 



tions for WIMP-nucleon interaction, a{x — n) [S], or 
WIMP annihilations [2] . Therefore the effective field the- 
ories allow us to map measurements performed at the 
LHC to the quantities relevant for direct-detection and 
indirect-detection dark matter search experiments. 

The effective-field-theory approach is valid as long as 
the unknown new mediator particles that couple the 
dark-matter particles to SM quarks or gluons are too 
heavy to be resolved: q < M* , where q is the momentum 
transfer. The breakdown of the effective approach de- 
pends ultimately on the details of the new and unknown 
physics, specifically on the number of new mediator par- 
ticles and the new couplings. Therefore, these theories 
cannot be treated generically and must be interpretted 
with some care. To guide the interpretation, we indicate 
the range of validity as lower bounds on the mass suppres- 
sion scale Mi, following ref. [3]. Assuming the simplest 
possible structure of new physics (mediation via exactly 
one new heavy mediator of mass M, M^, = Mj ^g\gi-, 
gi and g2 being coupling constants), bounds on the sup- 
pression scale can be placed by requiring M > 2m^ and 
that the new physics be as strongly coupled as possible 
for it to be still perturbative {^Jgig2 < 47r): 

Tn 

Mi, > -^(D5 to D14 and C3 to C6), 
27r 

\ —>^ pltoD4 , 
\ m„ 27r 



^ > ^ (CI and C2). 



Note that we are accounting for additional factors of rUq 
in the definitions of operators Dl to D4 and CI, C2 of 
ref. 



Experimental Searches 

The experimental searches typically require one or 
more high-px object and missing transverse momentum, 
see Table |T] for a summary and comparison of the mono- 
photon and mono-jet selections. 

The mono-Z analysis [TU] uses the ATLAS ZZ livv 
cross-section measurement [5], which requires: 

• two same-flavor opposite-sign electrons or muons, 
each with p^r^ > 20 GeV, \ri^\ < 2.5; 

• dilepton invariant mass close to the Z boson mass: 
ma G [mz — 15, mz + 15] GeV; 

• no particle-level jet with > 25 GeV and 
W\ <4.5; 



(bf-p|l)/p| <0.6; 

-p!^"" X cos(A0(p^'',pf )) > 80 GeV. 
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TABLE I; Summary of event selection requirements in 
ATLAS and CMS mono-jet or mono-photon analyses. Note 
that ATLAS uses two signal regions {j^t > 350 or 500 GeV) 
for the mono-jet analyses, depending on the operator. 





ATLAS 


CMS 


jet 


1 or 2 jets 


1 or 2 jets 




pi^ > 350(500) GeV 


> 110 GeV 




p;^^ > 30 GeV 


p!^ > 30 GeV 




PT > 350(500J GeV 


pT > 350 GeV 




veto leptons 


veto leptons 




A<j}{j2,^T) > 0.5 


A0(ji,j2) <2.5 


7 


1 photon, pt > 150 GeV 


1 photon Pt > 145 GeV 




> 150 GeV 


> 130 GeV 




< 1 jet with Pt > 30 GeV 


track with pr > 20 GeV 




isolation details 


isolation details 




A<?!)(7,^t) > 0.4 
A(^(ji,^T) > 0.4 






veto leptons 





no The selection efficiency of each selection for each opera- 

111 tor is given in Table [TI| and were estimated in the following 

112 way. References [1H7] provide signal efficiency for several 

113 select operators; this efficiency is the product of geomet- 

114 ric and kinematic acceptance of the selection criteria and 

115 object reconstruction efficiency. The object reconstruc- 

116 tion efficiency depends on the details of the detector per- 
il? formance, but is largely independent of operator. The 
lis geometric and kinematic acceptances can be reliably es- 

119 timated using parton-level simulated event samples. We 

120 measure the geometric and kinematic efficiency for each 

121 operator, and use the quoted total efhciences to deduce 

122 the object reconstruction efficiencies. This allows us to 

123 estimate the total efficiency for each operator. 



Combination 

The separate analyses, each of which are single-bin 
counting experiments, are combined into a multi-bin 
counting experiment. This allows for a coherent signal 
rate to be tested across channels, but preserves their dis- 
tinct signal-to-background ratios. 

The background estimates are taken directly from the 
experimental publications, see a summary in Table |III[ 
and are assumed to be uncorrelated across channels, 
as they are typically dominated by channel-specific or 
detector-specific uncertainties. For example, in some 
cases the background estimates are data-driven, and the 
dominant uncertainties are in the finite statistics of in- 
dependent control samples. Inclusion of correlations up 
to 20% does not qualitatively impact the results of the 
combination. 

The backgrounds, their uncertainties and the observed 
yield can be used to calculate a 90% CL upper limit on' 
the number of signal events N in the sample, see TablelTTT] 



TABLE II: Selection efficiency as percentages for each channel 
of the analyses used in the combination, for operators Dl- 
14 and C1-C6 for low and high values of the WIMP mass 
rrix- ATLAS mono-jet analysis has two signal regions, we 
use J^T > 500(350) GeV and pi^ > 500(350) GeV region 
for operators D9-D14 (D1-D8 and C1-C6). Operators Dll- 
14, C5 and C6 only couple to gluon initial states, and so 
have no efficiency for photon or Z boson radiation. The Z 
efficiencies include the Z U branching fraction. Jet and 
photon samples include a pt > 80 GeV cut at parton-level. 







ATLAS 


CMS 


Operator 


jet 


7 


Z 


jet 


7 


Dl 


10 


0.4% 


11.2% 


1.2% 


0.7% 


8.0% 




1000 


2.6% 


19.1% 


1.2% 


3.6% 


11.3% 


D2 


10 


0.4% 


10.8% 


1.2% 


0.7% 


8.0% 




1000 


2.4% 


18.6% 


1.1% 


3.7% 


11.3% 


D3 


10 


0.5% 


11.1% 


1.2% 


0.7% 


8.0% 




1000 


2.6% 


18.9% 


1.2% 


3.9% 


11.3% 


D4 


10 


0.5% 


10.8% 


1.2% 


0.7% 


7.6% 




1000 


2.6% 


18.6% 


1.1% 


3.7% 


11.3% 


D5 


10 


1.7% 


18.2% 


0.9% 


2.2% 


11.3% 




1000 


3.3% 


23.5% 


1.1% 


4.5% 


14.7% 


D6 


10 


1.7% 


18.7% 


0.9% 


2.2% 


12.0% 




1000 


3.2% 


23.6% 


1.1% 


4.4% 


15.2% 


D7 


10 


1.7% 


18.1% 


0.9% 


2.4% 


11.3% 




1000 


3.3% 


23.4% 


1.1% 


4.4% 


14.5% 


D8 


10 


1.7% 


18.5% 


0.9% 


2.3% 


11.8% 




1000 


3.1% 


23.6% 


1.1% 


4.3% 


15.1% 


D9 


10 


0.9% 


23.5% 


1.4% 


4.1% 


14.1% 




1000 


1.2% 


23.3% 


1.4% 


5.1% 


14.8% 


DIO 


10 


1.1% 


23.6% 


1.4% 


4.2% 


14.4% 




1000 


1.2% 


23.4% 


1.4% 


5.2% 


14.8% 


Dll 


10 


0.9% 


- 


- 


4.1% 


- 




1000 


2.4% 


— 


— 


7.5% 


— 


D12 


10 


1.0% 


— 


— 


4.2% 


— 




1000 








i .4/0 




D13 


10 


0.9% 






4.1% 






1000 


2.4% 






7.5% 




D14 


10 


1.1% 






4.0% 






1000 


2.4% 






7.4% 




CI 


10 


0.1% 


7.0% 


1.0% 


0.2% 


5.3% 




1000 


2.3% 


18.2% 


1.1% 


3.3% 


11.0% 


C2 


10 


0.1% 


7.0% 


1.0% 


0.1% 


5.6% 




1000 


2.5% 


18.4% 


1.1% 


3.8% 


11.2% 


C3 


10 


1.7% 


18.4% 


0.9% 


2.3% 


11.6% 




1000 


2.9% 


23.6% 


1.1% 


4.1% 


14.9% 


C4 


10 


1.4% 


18.4% 


0.9% 


2.2% 


11.8% 




1000 


3.0% 


23.8% 


1.1% 


4.1% 


15.3% 


05 


10 


1.4% 






1.7% 






1000 


5.9% 






7.6% 




06 


10 


1.2% 






1.7% 






1000 


5.9% 






7.6% 





is almost completely model independent. Translating it 
into a limit on the cross section for the pp — )■ xx + ^ sig- 
nal requires the effieciency of the signal in each selection. 



and Table IV using the CLs method [Ulllj. This valueus 



see Table III These individual limits reproduce well the 
results reported by the experiments. 

The signal regions are nearly orthogonal, but not ex- 
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TABLE III: 90% CL limits on iVovonts, efficiencies for = 10 
GeV, and limits on cripp — >■ XX + using the D5 operator. 
In the case of the Z + final state, the efficiency is relative 
to Z ^ ££ decays only. 



TABLE V: 90% CL hmits on a{pp XX + X) for = 10 
GeV, theory prediction for A/* — 1 TeV, and limits on M* 
using the D5 operator. In the case of the Z + ]^t final state, 
the prections include the Z ^ U branching fraction. 



Channel 


Bg. 


Obs 


Limit Eff. 


Lumi. 


Limit 


Channel 


Limit a 


Pred. Limit Af* 










N 


(tb-i) 


a (fb) 




(fb) 


(fb) 


(GeV) 




ATLAS jet+^T 


750 ± 60 


785 


139.3 1.7% 


4.8 


1,700 


ATLAS jet±_^T 


1,700 


370 


685 -1 

750 J ''S^ -| 




CMS jet+^T 


1225 ± 101 


1142 


125.2 2.2% 


5.0 


1,140 


CMS jet±^T 


1,140 


370 




ATLAS 7 + 


137 ± 20 


116 


27.4 18% 


4.6 


33 


ATLAS 7 ± 


33 


3.7 


580 T 


^ 795 


CMS 7 + 


75.1 ±9.4 


73 


19.3 11% 


5.0 


35 


CMS 7±^T 


35 


3.7 


570 1 645 J 




ATLAS Z + 


86.2 ±7.2 


87 


21.7 13% 


4.6 


36 


ATLAS Z ± 


36 


0.5 


340 





TABLE IV: 90% CL limits on iVcvents, efficiencies for = 10. 
GeV, and limits on (t(pp — XX ± X) using the D9 operator. 



Channel 


Bg. 


Obs 


Limit Eff. Lumi. 


Limit 








TV (fb-i) 


a (fb) 


ATLAS jet±| 


?T 83 ± 14 


77 


25.5 0.9% 4.8 


590 


CMS jet±^T 


1225 ± 101 


1142 


125.2 4.1% 5.0 


610 



190 



150 actly. For example, the mono-jet analyses do not veto 

192 

151 events with a photon, and the mono-photon analyses al- 

193 

152 low the presence of one jet. From our parton-level sim- 

153 ulated event samples, we estimated the overlaps among 

195 

154 different channels and found that the overlap fraction is 

155 less than 1%. 

156 The individual analyses include signal uncertainties of 

(W . .196 

157 up to 20% on the cross section, mostly due to uncertain- 

158 ties in jet energy calibration and levels of initial-state 

159 radiation. These uncertainties do not affect the cross-"^ 

160 section limits, but can be simply applied to limits on™ 

161 Mi^. In each case, we quote the limit using the central™ 

162 value. 

163 To summarize, the assumptions made in this combina-^"^ 

164 tion are 

203 

165 • the background uncertainties are monolithic and204 

166 uncorrelated, and 205 

206 

167 • the signal selections are orthogonal 

168 Combining channels is then straightforward, though 

169 the intermediate step of a model-independent limit on^°^ 

170 the number of events N is no longer possible, as the lim- 

171 its depend on the relative distribution of signal eventS208 

172 across channels, which is model specific. Instead, cross-209 

173 section limits are obtained directly. These limits are2io 

174 then converted into limits on M^, using the relationships2ii 

175 from Ref. [14j . The individual-channel limits, combina-212 

176 tion across experiments and the grand combination of alLu 

177 channels are shown in Table |V] for the D5 operator and2i4 

178 one choice of m^. Clearly the mono-jet analyses are the2i5 

179 most powerful, and the greatest gain in combination is2i6 

180 from combining the ATLAS and CMS mono-jet analy-217 

181 ses, though the addition of the mono-photon and mono-218 

182 Z gives a non-negligible improvement in the combined2i9 

183 result. 220 



Limits on for the D5 and D8 operators are shown in 
Fig[3]and|4] as well as limits on cr(x — Where the M^, 
limits exceed the thermal relic values taken from Ref. [3J , 
assuming that dark matter is entirely composed of ther- 
mal relics, the resulting dark matter density of the uni- 
verse would contradict WMAP measurements; therefore, 
WIMPs cannot couple to quarks or gluons exclusively via 
the given operator and account entirely for the relic den- 
sity. This my. region is either excluded, or requires that 
annihilation channels to leptons must exist, or partici- 
pation of different operators which interfere negatively, 
thereby reducing the limits on Mi,. 

Application to other models 

While the experimental results are usually quoted for 
a small selection of the effective operator models, the 
analyses are clearly relevant for all of them. 

We re-interpret the experimental analyses in the con- 
text of each operator and perform the grand combination 
across all channels. Figure [5] and Table [Vl] show the lim- 
its on Mi,, translated to the WIMP-nucleon cross section 
where possible. In addition, we translate the limits on 
D5 and D8 into limits on the WIMP annihilation cross 
section, see Fig. [6] 

Conclusions 

We have presented the first combination of collider- 
based searches for dark matter pair production, us- 
ing final states involving jets, photons and leptonically- 
decaying Z bosons in the context of effective field the- 
ories. The most powerful results are from the mono-jet 
analyses, and the greatest gains come from the combina- 
tion of the independent analyses from ATLAS and CMS, 
though the other final states make a non-negligblc im- 
provement. The results are the strongest limits to date 
from collider searches in the effective field theory context. 

In addition, we have reinterpreted the experimental re- 
sults, quoted by ATLAS and CMS only for a few effective 
operators, across a broad range of operators, providing 
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FIG. 3: Limits at 90% CL in M* (top) and in the spin- 
independent WIMP-nucleon cross section (bottom) for indi- 
vidual and combined limits using the D5 operator as a func- 
tion of m^. 



FIG. 4: Limits at 90% CL in (top) and in the spin- 
dependent WIMP-nucleon cross section (bottom) for individ- 
ual and combined limits using the D8 operator as a function 
of m,),. 



221 a comprehensive view of the power of these searches to 

222 constrain the weak-level or weaker interactions between 

223 dark matter and standard model particles. 

236 

224 We have made use of the effective field theory frame-^j^ 

225 work to convert the ATLAS and CMS results to quan-238 

226 titles relevant for direct-detection and indirect-detection239 

227 dark matter searches. Under the assumptions made for^io 

228 the effective operators, LHC limits can be very compet-^''^ 

229 itive, in particular for low-mass dark matter particles^"^ 

230 < 10 GeV. ^"^ 
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FIG. 5: Combined limits on M* at 90% CL, using all avail- 
able channels, for operators Dl-14 and C1-C5 as a function 
of my. 
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FIG. 6: Top and center, limits at 90% CL on the spin- 
independent and spin-dependent WIMP-nucleon cross sec- 
tion, — n.), for available operators. Bottom, interpretation 
of the limits on D5 and D8 in terms of the velocity-averaged 
WIMP-annhiliation cross section, as defined in Ref [2]. 
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TABLE VI: Combined limits on at 90% CL, using all 
available channels, for operators Dl-11 and C1-C5 for low and 
high values of the WIMP mass m^. Where possible, limits 
are shown on the WIMP-nucleon cross section, a{x — n). 



Operator 


m;^ 


M* 


<^{x - 


n) 




(GeV) (GeV) 


[cm^] 


Dl 


10 


34 


5.2x10" 


-ay 




1000 


10 


8.3x10" 


-36 


D2 


10 
1000 


34 
13 






D3 


10 
1000 


34 
10 






D4 


10 
1000 


34 
13 






D5 


10 


795 


2.0x10" 


-39 




1000 


325 


8.8x10" 


-38 


D6 


10 
1000 


791 
221 






D7 


10 
1000 


812 
324 






D8 


10 


811 


6.5x10" 


-4i 




1000 


222 


1.4x10" 


-38 


D9 


10 


1331 


8.9x10" 


-42 




1000 


413 


1.1x10" 


-39 


DIO 


10 
1000 


1410 
415 






Dll 


10 


339 


1.8x10" 


-44 




1000 


155 


2.4x10" 


-42 


D12 


10 
1000 


342 
188 






D13 


10 
1000 


427 
195 






D14 


10 
1000 


429 
237 






CI 


10 


8 


4.2x10" 


-36 




1000 


1 


6.1x10" 


-36 


O 


10 
1000 


8 
1 






C3 


10 


575 


7.5x10" 


-39 




1000 


153 


1.8x10" 


-36 


C4 


10 
1000 


556 
154 






C5 


10 


201 


4.4x10" 


-4i 




1000 


41 


3.1x10" 


-42 


C6 


10 
1000 


286 
57 
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Appendix: Individual Operators 



285 In Figs. [7) m [9) [TO][TTJ and [12] we show the combined^''^ 

286 hmits for each operator, compared to the thermal rehc^'^ 

287 values. Where the limits exceed the thermal relic val-^'" 

288 ues, assuming that dark matter is entirely composed oP'^ 

289 thermal relics, the dark matter density of the universe 



would contradict measurements and hence cannot couple 
to quarks or gluons exclusively via the given operator. 
This region is either excluded, or else other annihi- 
lation channels to leptons must exist, or finally different 
operators may interfere negatively thereby reducing the 
limits on M-^. 



m, [GaV] 








FIG. 7: Combined limits on versus dark matter mass 
for operators Dl, D2, D3 and D4. The Mi, values at which 
dark matter particles of a given mass would result in the re- 
quired relic abundance are shown as green lines [3], assuming 
annihilation in the early universe proceeded exclusively via 
the given operator. 



m, [GaV] 







/ - 








10 10^ 10 

[GeV] 


/ 













[GaV] 



FIG. 8: Combined limits on M^, versus dark matter mass 
for operators D5, D6, and D7. The M* values at which 
dark matter particles of a given mass would result in the re- 
quired relic abundance are shown as green lines [3], assuming 
annihilation in the early universe proceeded exclusively via 
the given operator. 




FIG. 9: Combined limits on M^, versus dark matter mass 
for operators D8,D9 and DIO. The M* values at which 
dark matter particles of a given mass would result in the re- 
quired relic abundance are shown as green lines [3], assuming 
annihilation in the early universe proceeded exclusively via 
the given operator. 
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m, [GeV] 



FIG. 10: Combined limits on versus dark matter mass 
for operators Dll, D12, D13 and D14. The M* values 
at which dark matter particles of a given mass would result 
in the required relic abundance are shown as green lines [3], 
assuming annihilation in the early universe proceeded exclu- 
sively via the given operator. 
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FIG. 11: Combined limits on M* versus dark matter mass 
for operators CI, C2, and C3. The Af* values at which 
dark matter particles of a given mass would result in the re- 
quired relic abundance are shown as green lines [3], assuming 
annihilation in the early universe proceeded exclusively via 
the given operator. 
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m, [GaV] 




m, [GeV] 



FIG. 12: Combined limits on M* versus dark matter mass 
for operators C4, C5 and C6. The M* values at which 
dark matter particles of a given mass would result in the re- 
quired relic abundance are shown as green lines [3], assuming 
annihilation in the early universe proceeded exclusively via 
the given operator. 



